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Abstract. First results of a numerical investigation of primordial black hole formation in the 
■ radiation dominated phase of the Early Universe are presented. The simulations follow the 
gravitational collapse of three different families of high-amplitude density fluctuations imposed 
at the time of horizon crossing. The threshold for black hole formation, <5 C ~ 0.7, is found to be 
\& ' nearly identical for all perturbation families if the control parameter, S, is chosen as the total 
excess mass within the initial horizon volume. Furthermore, we demonstrate that the scaling 
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of black hole mass with distance from the formation threshold, known to occur in near-critical 
| gravitational collapse, applies to primordial black hole formation. 

q 1. Introduction 

h ; 

Every quantitative analysis of the primordial black hole (PBH) number and mass spectrum [ 
H| requires knowledge of the threshold parameter, <5 C , separating perturbations that form black 
. holes from those that do not, and the resulting black hole mass, Mbh, as a function of distance 
from the threshold. In order to determine S c and Mbh for various initial conditions, we performed 
one-dimensional, general relativistic simulations of the hydrodynamics of PBH formation in the 
radiation-dominated phase of the early universe. Three families of perturbation shapes were 
chosen to represent "generic" classes of initial data, reflecting the lack of information about 
the specific shape of primordial fluctuations. The numerical technique is sketched in Section 
^. Defined as the excess mass within the horizon sphere at the onset of the collapse, we find 
5 C ~ 0.7 for all three perturbation shapes, indicating that the threshold value may indeed be 
universal (Section ||). A numerical confirmation of the previously suggested power-law scaling of 
Mbh with 5 — 5 C [ , related to the well-known behavior of collapsing space-times at the critical 
point of black hole formation [ |j| , is presented in Section |||. In this framework, the PBH mass 
spectrum is determined by the dimensionless coefficient, K, the scaling exponent, 7, and the 
initial horizon mass, Mh, such that 

M bh = KM h {6 - S c y . (1) 

We provide numerical results for K and 7 for the three perturbation families. A more detailed 
description of the numerical technique and further results will be presented in a forthcoming 
publication [ . 
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2. Numerical technique 

The dynamics of collapsing density perturbations in the Early Universe are determined by 
the general relativistic equations of motion for a perfect fluid, the field equations, the first law 
of thermodynamics, and a radiation-dominated equation of state. The assumption of spherical 
symmetry is well justified for large fluctuations in a Gaussian distribution [ ||, reducing the 
problem to one spatial dimension. 

For our simulations, we chose the formulation of the hydrodynamical equations by Hernandez 
and Misner [ |(| as implemented by Baumgarte et al. [ Q . Based on the original equations by Mis- 
ner and Sharp [ g], Hernandez and Misner proposed to exchange the Schwarzschild time variable, 
t, with the outgoing null coordinate, u. In so doing, the hydrodynamical equations retain the 
Lagrangian character of the Misner-Sharp equations but avoid crossing into the event horizon 
of a black hole once it has formed. Covering the entire space-time outside while asymptotically 
approaching the event horizon, the Hernandez-Misner equations are perfectly suited to follow 
the evolution of a black hole for long times after its formation without encountering coordinate 
singularities. This allowed us, in principle, to study the accretion of material onto newly formed 
PBHs for arbitrarily long times (in contrast with earlier calculations [ [|). Since the expanding 
outer regions of our simulated piece of the universe are most conveniently tracked in a comoving 
numerical reference frame, the Lagrangian form of the Hernandez-Misner equations is their sec- 
ond major asset. It also provides a simple prescription for the outer boundary condition, which 
is defined to match the exact solution of the Friedmann equations for a radiation dominated flat 
universe. Hence, the pressure follows the analytic solution 

W 

where r is the proper time of the outermost fluid element (corresponding to the cosmological 
time t in a Friedmann-Robertson- Walker (FRW) universe) and Po an d To are the initial values 
for pressure and proper time. 

3. Threshold for black hole formation 

We studied the spherically symmetric evolution of three families of curvature perturbations. 
Initial conditions were chosen to be perturbations of the energy density, e = poe, in unperturbed 
Hubble flow specified at horizon crossing. The first family of perturbations is described by a 
Gaussian-shaped overdensity that asymptotically approaches the FRW solution at large radii. 
The other two families of initial conditions involve a spherical Mexican Hat function and a fourth 
order polynomial. These functions are characterized by rarefaction regions outside of the horizon 
radius, i?h, that identically compensate for the additional mass of the overdensities inside the 
horizon volume, so that the mass derived from the total integrated density profile is equal to 
that of an unperturbed FRW solution. In our numerical experiments, the amplitude, A, of the 
perturbations is used to tune the initial conditions to sub- or supercriticality with respect to 
black hole formation. The critical amplitude, A c , is strongly shape-dependent, varying between 
A c = 3.04 for Mexican-Hat-shaped perturbations and A c = 2.05 for the Gaussian curve. If, 
however, we define the control parameter 5 as the additional mass inside i?h in units of the 
horizon mass, we find strikingly similar values — 5 C = 0.67 (Mexican Hat), 5 C = 0.70 (Gaussian 
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curve), and 5 C = 0.71 (polynomial) — for all three families of initial data in our study. This 
number is considerably greater than the previously employed threshold <5 C = 1/3 following from 
analytical estimates [ [l| . Given the qualitative difference of the functional forms of the different 
perturbation families, the result that 5 C of the Gaussian perturbation lies in between the critical 
values of the mass compensated functions is surprising. Our results suggest that 5 C ~ 0.7 
(with 5 defined as above) is a universal statement, i.e., true for all perturbation shapes in the 
radiation-dominated regime. This remains to be verified by means of additional experiments. 



4. Scaling of PBH masses with distance to the threshold 

For a variety of matter models, it is well-known that the dynamics of near-critical collapse 
exhibit continuous or discrete self-similarity and power-law scaling of the black hole mass with 
the offset from the critical point {Eq. (jlj) [|3|, [l0|]}. In particular, Evans and Coleman [ 11] found 



self-similarity and mass scaling in numerical experiments of a collapsing radiation fluid. They 
numerically determined the scaling exponent 7 ~ 0.36, followed by a linear perturbation analysis 



of the critical solution by Koike et al. [ 12 1 that yielded 7 ~ 0.3558. Until recently, it was believed 
that entering the scaling regime requires a degree of fine-tuning of the initial data that is un- 
natural for any astrophysical application. It was noted [ [2j that fine-tuning to criticality occurs 
naturally in the case of PBHs forming from a steeply declining distribution of primordial den- 
sity fluctuations, as generically predicted by inflationary scenarios. In the radiation-dominated 
cosmological epoch, the only difference with the fluid collapse studied by Evans and Coleman [ 



11 1 is the expanding, finite-density-background space-time of a FRW universe. Assuming that 
self-similarity and mass scaling are consequences of an intermediate asymptotic solution that 
is independent of the asymptotic boundary conditions, Eq. (|l|) is applicable to PBH masses, 
allowing the derivation of a universal PBH initial mass function [ |2| . Figure (|l|) presents numer- 
ical evidence for mass scaling according to Eq. (||) in black hole collapse in an asymptotic FRW 
space-time. All three perturbation families give rise to scaling solutions with a scaling exponent 
7 w 0.36. 



5. Conclusions 

This work discusses numerical collapse simulations of three generic families of energy density 
perturbations, one with a finite total excess mass with respect to the unperturbed FRW solution 
and two mass compensated ones. Among various possible definitions for the collapse control 
parameter 6, the total excess gravitational mass of the perturbed space-time with respect to the 
unperturbed FRW background enclosed in the initial horizon volume is the only one that gives 
rise to a similar threshold value for all three shape families, 5 C ~ 0.7. Whether this result is an 
indication for universality of 5 C in this specific definition needs to be verified with the help of 
additional simulations using a larger sample of initial perturbation shapes. 

The previously suggested [ |2| scaling relation between Mbh and 5 — 5 C , based on the analogy 
with critical phenomena observed in near-critical black hole collapse in asymptotically non- 
expanding space-times [ ||, is confirmed numerically for an asymptotic FRW background. For 
the smallest black holes in this investigation, the scaling exponent is 7 ~ 0.36, which is identical 



to the non-expanding numerical and analytical results [ 11, 12 1 within our numerical accuracy. 
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Figure 1. Black hole masses as a function of 5— 5 C for three different perturbation shape families. 
The best fit parameters to equation (Q) are: 7 = 0.36, K = 2.85, 5 C = 0.6745 (Mexican Hat 
perturbation, triangles); 7 = 0.37, K = 2.39, 5 C = 0.7122 (polynomial perturbation, crosses); 
7 = 0.34, K = 11.9, 5 C = 0.7015 (Gaussian curve perturbation, diamonds). 



The parameter K of Eq. (|l]), needed in addition to 7 to evaluate the two-parameter PBH IMF 
derived in [ |2j, ranges from K ~ 2.4 to K ~ 12. 

The author wishes to thank K. Jedamzik for a fruitful collaboration, T. Baumgarte for pro- 
viding the original version of the hydrodynamical code, and Joan George for valuable stylistic 
corrections. 
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